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Abstract 



We study random graphs with an i.i.d. degree sequence of which the tail of the distribution 
function F is regularly varying with exponent r € (1, 2). Thus, the degrees have infinite mean. 
Such random graphs can serve as models for complex networks where degree power laws are 
[V] observed. 

P I ■ The minimal number of edges between two arbitrary nodes, also called the graph distance 

or the hopcount, in a graph with N nodes is investigated when N — > oo. The paper is part of 
a sequel of three papers. The other two papers study the case where r 6 (2, 3), and r € (3, oo), 
respectively. 

The main result of this paper is that the graph distance converges for t £ (1,2) to a limit 
random variable with probability mass exclusively on the points 2 and 3. We also consider the 
case where we condition the degrees to be at most N a for some a > 0. For r _1 < a < (r — 
the hopcount converges to 3 in probability, while for a > (t — 1) , the hopcount converges to 
the same limit as for the unconditioned degrees. Our results give convincing asymptotics for 
\ the hopcount when the mean degree is infinite, using extreme value theory. 

O . Key words and phrases: Extreme value theorem, Internet, random graphs. 
q . AMS 2000 Subject classifications: Primary 60G70, Secondary 05C80. 

1 Introduction 

The study of complex networks has attracted considerable attention in the past decade. There are 
numerous examples of complex networks, such as social relations, the World-Wide Web and Inter- 
net, co-authorship and citation networks of scientists, etc. The topological structure of networks 
affects the performance in those networks. For instance, the topology of social networks affects 
the spread of information and disease (see e.g., |19j). while the performance of traffic in Internet 
depends heavily on the topology of the Internet. 

Measurements on complex networks have shown that many real networks have similar prop- 
erties. A first example of such a fundamental network property is the fact that typical distances 
between nodes are small. This is called the 'small world' phenomenon, see the pioneering work of 
Watts |20j and the references therein. In Internet, for example, e-mail messages cannot use more 
than a threshold of physical links, and if the distances in Internet would be large, e-mail service 
would simply break down. Thus, the graph of the Internet has evolved in such a way that typical 
distances are relatively small, even though the Internet is rather large. 
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A second, maybe more surprising, property of many networks is that the number of nodes with 
degree k falls of as an inverse power of k. This is called a 'power law degree sequence'. In Internet, 
the power law degree sequence was first observed in j^j. The observation that many real networks 
have the above properties have incited a burst of activity in network modelling. Most of the models 
use random graphs as a way to model the uncertainty and the lack of regularity in real networks. 
See an d the references therein for an introduction to complex networks and many examples 
where the above two properties hold. 

The current paper presents a rigorous derivation for the random fluctuations of the distance 
between two arbitrary nodes (also called the geodesic) in a graph with i.i.d. degrees with infinite 
mean. The model with i.i.d. degrees is a variation of the configuration model, which was originally 
proposed by Newman, Strogatz and Watts where the degrees originate from a given deter- 
ministic sequence. The observed power exponents are in the range from r = 1.5 to r = 3.2 (see 
[HI Table II] or Table II]). In a previous paper with Van Mieghem we have investigated 
the case r > 3. The case r G (2, 3) was studied in [T5|. Here we focus on the case r G (1, 2), and 
study the typical distances between arbitrary connected nodes. In a forthcoming paper JH], we 
will survey the results from the different cases for r, and investigate the connected components of 
the random graphs. 

This section is organized as follows. In Section ll.ll we start by introducing the model, and in 
Section [1.21 we state our main results. In Section fl .31 we explain heuristically the results obtained. 
Finally we describe related work in Section fl .41 

1.1 The model 

Consider an i.i.d. sequence D\, D2, ■ ■ ■ , D N . Assume that L N = YljLi Dj is even. When L N is odd, 
then we increase the number of stubs of Djy by 1, i.e., we replace Dn by Djv + 1. This change will 
make hardly any difference in what follows, and we will ignore it in the sequel. 

We will construct a graph in which node j has degree Dj for all 1 < j < N. We will later 
specify the distribution of Dj. We start with N separate nodes and incident to node j, we have Dj 
stubs which still need to be connected to build the graph. 

The stubs are numbered in an arbitrary order from 1 to L N . We now continue by matching at 
random the first stub with one of the L N — 1 remaining stubs. Once paired, two stubs form an edge 
of the graph. Hence, a stub can be seen as the left or the right half of an edge. We continue the 
procedure of randomly choosing and pairing the next stub and so on, until all stubs are connected. 

The probability mass function and the distribution function, of the nodal degree, are denoted 

by 

F(D 1= j) = fj, j = 0,1,2,..., and F(z) = (1.1) 

j=0 

where |_^J is the largest integer smaller than or equal to x. Our main assumption will be that 

1 -F(x) =x- T+1 L(x), (1.2) 

where r G (1,2) and L is slowly varying at infinity. This means that the random variables Dj obey 
a power law with infinite mean. The factor L is meant to generalize the model. 

1.2 Main results 

We define the graph distance H N between the nodes 1 and 2 as the minimum number of edges 
that form a path from 1 to 2, where, by convention, this distance equals 00 if 1 and 2 are not 
connected. Observe that the distance between two randomly chosen nodes is equal in distribution 
to H N , because the nodes are exchangeable. 

We will present in this paper two separate theorems for the case r G (1,2) . In the first theorem 
we take the sequence D\,D2, . . . , D N an i.i.d. sequence with distribution F, satisfying ()1.2j) . with 
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r G (1,2). The result is that the graph distance or hopcount converges in distribution to a limit 
random variable with mass p = p T , 1 — p, on the values 2, 3, respectively. 

One might argue that including degrees larger than iV — 1 is artificial in a network with N 
nodes. In fact, in many real networks, the degree is bounded by a physical constant. Therefore 
we also consider the case where the degrees are conditioned to be smaller than N a , where a is an 
arbitrary positive number. Of course, we cannot condition on the degrees to be at most M, where 
M is fixed and independent on N, since in this case, the degrees are uniformly bounded, and this 
case is treated in Therefore, we consider cases where the degrees are conditioned to be at 

most a given power of N. 

The result with conditioned degrees appears in the second theorem below. It turns out that for 
a > 1/(t — 1), the conditioning has no influence in the sense that the limit random variable is the 
same as that for the unconditioned case. This is not so strange, since the maximal degree is of order 
jV 1 /( r - 1 ) ) so that the conditioning does nothing in this case. However, for 1/r < a < l/(r — 1), 
the graph distance converges to a degenerate limit random variable with mass 1 on the value 3. It 
would be of interest to extend the possible conditioning schemes, but we will not elaborate further 
on it here. 

Theorem 1.1 Fixr G (1,2) in M.Ujl and let D\, D2, ■ ■ ■ ,D N denote an i.i.d. sequence with common 
distribution F . Then, 

lim F(H N = 2) = 1 - lim F(H N = 3) = p, (1.3) 

where p = p T G (0, 1). 

In the theorem below we will write D ( - N) for the random variable D conditioned on D < N a . 
Thus, 

F(D (N) = k) = — — — r, 0<k<N a . (1.4) 

v ' F(D < N a ) ~ y ' 

Theorem 1.2 Fix r G (1, 2) in ll.2\) and let D[ ' ,D% , ■■■ , D^' be a sequence of i.i.d. copies of 
D (N) . 

(i) If 1/r < a < 1/(t - 1), then 

lim F(H N = 3) = 1. (1.5) 

N— >oo 

(ii) If a > l/(r - 1), then 

lim F(H N = 2) = 1 - lim F(H N = 3) = p, (1.6) 

N^oo N-too 

where p = p T G (0, 1). 

Remark: For a < 1/t, we have reasons to believe that F(H N > 3) remains uniformly positive 
when N — > 00. A heuristic for this fact is given in the next section. 

1.3 Heuristics 

When r G (1,2), we consider two different cases. In Theorem ll.il the degrees are not conditioned, 
while in Theorem 11.21 we condition on each node having a degree smaller than N a . We will now 
give a heuristic explanation of our results. 

In two previous papers we have treated the cases r G (2, 3) and r > 3. In these cases, 

it follows that the probability mass function {fj} introduced in (|l.lj) has a finite mean, and the 
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the number of nodes on graph distance k from node 1 can be coupled to the number of individuals 
in the k th generation of a branching process with offspring distribution {gj} given by 

9i = (1-7) 

where [i = E[Lq]. For r G (1,2), as we are currently investigating, we have \i = do, and the 
branching process used in |11[ I12| does not exist. 

When we do not condition on Dj being smaller than N a , then L N is the i.i.d. sum of N random 
variables D\,D2, ■ ■ ■ ,D N , with infinite mean. It is well known that in this case the bulk of the 
contribution to L N ~ JV 1 ' ( T— ^ comes from a finite number of nodes which have giant degrees (the 
so-called giant nodes). A basic fact in the configuration model is that two sets of stubs of sizes n 
and m are connected with high probability when nm is at least of order L N . Since the giant nodes 
have degree roughly N 1 ^ T ~ 1 > , which is much larger than y/L N , they are all attached to each other, 
thus forming a complete graph of giant nodes. Each stub is with probability close to 1 attached to 
a giant node, and therefore, the distance between any two nodes is, with large probability, at most 
3. In fact this distance equals 2 when the two nodes are connected to the same giant node, and is 
3 otherwise. 

When we truncate the distribution as in p.4|) . with a > 1/(t — 1), we hardly change anything 
since without truncation with probability 1 — o(l) all degrees are below N a . On the other hand, 
if a < 1/(t — 1) when, with truncation, the largest nodes have degree of order N a , and L N ~ 
jyl+a(2-r) _ Again, the bulk of the total degree L N comes from nodes with degree of the order 
N a , so that now these are the giant nodes. Hence, for 1/r < a < l/(r — 1), the largest nodes 
have degrees much larger than \/L N , and thus, with probability 1 — o(l), still constitute a complete 
graph. The number of giant nodes converges to infinity, as N — > oo. Therefore, the probability that 
two arbitrary nodes are connected to the same giant node converges to 0. Therefore, the hopcount 
equals 3 with probability converging to 1. If a < 1/r, then the giant nodes no longer constitute a 
complete graph suggesting that the resulting hopcount can be greater than 3. It remains unclear 
to us whether the hopcount converges to a single value (as in Theorem II. 2|) . or to more than one 
possible value (as in Theorem II. We do expect that the hopcount remains uniformly bounded. 

The proof in this paper is based on detailed asymptotics of the sum of N i.i.d. random variables 
with infinite mean, as well as on the scaling of the order statistics of such random variables. The 
scaling of these order statistics is crucial in the definition of the giant nodes which are described 
above. The above considerations are the basic idea in the proof of Theorem 11.11 In the proof of 
Theorem 11.21 we need to investigate what the conditioning does to the scaling of both the total 
degree L N , as well as to the largest degrees. 

1.4 Related work 

The above model is a variation of the configuration model. In the usual configuration model one 
often starts from a given deterministic degree sequence. In our model, the degree sequence is an 
i.i.d. sequence D±, . . . ,D N with distribution equal to a power law. The reason for this choice is 
that we are interested in models for which all nodes are exchangeable, and this is not the case when 
the degrees are fixed. The study of this variation of the configuration model was started in for 
the case r > 3 and studied by Norros and Reittu ^Hj i n case r S (2, 3). 

For a complete survey to complex networks, power law degree sequences and random graph 
models for such networks, see (3] and ^]. There a heuristic is given why the hopcount scales 
proportionally to logiV, which is originally from The argument uses a variation of the power 
law degree model, namely, a model where an exponential cut off is present. An example of such a 
degree distribution is 

fk = Ck- T e- k / K (1.8) 
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for some k > 0. The size of k indicates up to what degree the power law still holds, and where 
the exponential cut off starts to set in. The above model is treated in for any k < oo, but, for 
k = oo, falls within the regimes where r G (2,3) in |12j and within the regime in this paper for 
r G (1, 2). In [Uj, the authors conclude that since the limit as k — > oo does not seem to converge, 
the 'average distance is not well-defined when k < 3'. In this paper, as well as in ^2]> we show 
that the average distance is well-defined, but it scales differently from the case where r > 3. 

In the paper |13| . we give a survey to the results for the hopcount in the three different regimes 
r G (1,2), r G (2,3) and r > 3. There, we also prove results for the connectivity properties of 
the random graph in these cases. These results assume that the expected degree is larger than 2. 
This is always the case when r G (1, 2), and stronger results have been shown there. We prove that 
the largest connected component has size N(l + o(l)) with probability converging to one. When 
r G (1, |), we can even prove that with large probability, the graph is connected. When r > |, 
this is not true, and we investigate the structure of the remaining 'dust' that does not belong to 
the largest connected component. In the analysis, we will make use of the results obtained in this 
paper for r G (1,2). For instance, it will be crucial that the probability that two arbitrary nodes 
are connected converges to 1. 

There is substantial related work on the configuration model for the case r G (2,3) and r > 3. 
References are included in the paper ^2] for the case r G (2, 3) and in |TJ for r > 3. We again 
refer to the references in ^3] and f° r more details. The graph distance for r G (1,2), that 

we study here, has to our best knowledge not been studied before. Values of r G (1,2) have been 
observed in networks of e-mail messages and networks where the nodes consist of software packages 
(see ^3 Table II]), for which our configuration model with r G (1,2) can possibly give a good 
model. 

In pQ, random graphs are considered with a degree sequence that is precisely equal to a power 
law, meaning that the number of nodes with degree k is precisely proportional to k~ T . Aiello et 
al. pQ show that the largest connected component is of the order of the size of the graph when 
r < ro = 3.47875 . . ., where tq is the solution of Q{t — 2) — 2((t — 1) = 0, and where C is the 
Riemann Zeta function. When r > tq, the largest connected component is of smaller order than 
the size of the graph and more precise bounds are given for the largest connected component. When 
r G (1,2), the graph is with high probability connected. The proofs of these facts use couplings 
with branching processes and strengthen previous results due to Molloy and Reed • See also 

for a history of the problem and references predating |14l I15j . The problem of distances in the 
configuration model with r G (1,2) has, up to our best knowledge, not been addressed. See [2] for 
an introduction to the mathematical results of various models for complex networks (also called 
massive graphs), as well as a detailed account of the results in pQ. 

A detailed account for a related model can be found in jS] and [7], where links between nodes 
i and j are present with probability equal to WiWj / J2i w i fo r some 'expected degree vector' w = 
(wi, . . . ,W]\r). Chung and Lu [Hj show that when Wi is proportional to i~^r the average distance 
between pairs of nodes is log„ N(l + o(l)) when r > 3, and 2j]^rM^y(l + o(l)) when r G (2, 3). In 
their model, also r < 1 is possible, and in this case, similarly to r G (1, |) in our paper, the graph 
is connected with high probability. 

The difference between this model and ours is that the nodes are not exchangeable in j^j, but 
the observed phenomena are similar. This result can be understood as follows. Firstly, the actual 
degree vector in jE] should be close to the expected degree vector. Secondly, for the expected degree 
vector, we can compute that the number of components for which the degree is less than or equal 
to k equals 

\{i : Wi < k}\ cx \{i : i~^ < k}\ w £T T+1 . 

Thus, one expects that the number of nodes with degree at most k decreases as k~ T+1 , similarly 
as in our model. In 1 7] , Chung and Lu study the sizes of the connected components in the above 
model. The advantage of working with the 'expected degree model' is that different links are present 
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independently of each other, with makes this model closer to the random graph G(p,N). 
1.5 Organization of the paper 

The main body of the paper consists of the proofs of Theorem 11.11 in Section |2] and the proof of 
Theorem ll.2l in Section|3J Both proofs contain a technical lemma and in order to make the argument 
more transparent, we have postponed the proofs of these lemmas to the appendix. Section^Jcontains 
simulation results and some conclusions. 



2 Proof of Theorem 11.11 

In this section, we prove Theorem ll.il which states that the hopcount between two arbitrary nodes 
has, with probability 1 — o(l), a non-trivial distribution on 2 and 3. 

We will use an auxiliary lemma, which is a modification of the extreme value theorem for the 
k largest degrees, k G N. We introduce 

Ad <D { 2) <---<D (n) , 

to be the order statistics of D\, . . . , D N , so that D w = min{Di, . . . , D N }, D (2) is the second smallest 
degree, etc. Let (u N ) be an increasing sequence such that 

lim N(l- F(u N )) = 1. (2.1) 

N^oo 

It is well known that the order statistics of the degrees, as well as the total degree, are governed 
by u N in the case that r G (1,2). The following lemma shows this in some detail. 



Lemma 2.1 (a) for any k G N, 



'' ' 1 (£i, . . . , £fc) , in distribution, as N — > oo, (2-2) 



, . . . , 



where (£i, . . . , is a random vector with marginals in (0, oo) and with joint distribution function 
given, for any tuple < < ■ ■ ■ < y\ < oo, by 

IP (6 <yi,...,& <y k ) (2.3) 



E 



.l-T 



r\\ r 2 \ r k \ 

0<ri<-<ri-<fe 



-e 



-v k 



Moreover, 



£k — > 0, in probability, as k — > oo. (2-4) 



(b) 

where n is a random variable on (0, oo 



— — — > i], in distribution, as N — > oo, 
u N 



Proof. For part (a), we take pi = y\ T and u N {pi) = y%u N , i G {l,...,k}. Since 
L (N)N— (see e.g., 0) for some slowly varying function L (N), it follows from ()2.1j) that, 



lim N (1 - F(u N ( Pi ))) = Pl , i G {1, . . . , k}. 

N-+oo 
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Hence by 8, Theorem 4.2.6 and (4.2.4)], we have 

Pl)l ■ ■ ■ ! ^(JV-fc- 

pi 1 (p2 - piY 2 (pk - pk-iY k 



lim F(D m <u N (pi),...,D 

(N-k+l) < u N (pk)) 



E 

0<ri<-<r k <k 

y 

^— ' n! r?! rj.! 

0<ri<---<r fe <fc 



yt T)ri iy\- T - v\- T V ivl~ T ~ vl-lY" 



We compute now the marginal distribution of for k > 1. For any x > 0, due to (|2.3|) . we have 
P(£fc<x) =P(fi < oo,...,^ fc _i < oo,£ fc < x) 

k ( {l-r)n _ (l-r)r<_i\ 

= hm y e-^TT— — — - 

W>-,2/fc-l->oo ^ rj! 

0<ri<---<rfe_ 1 <rj.<fe i=l 



r( 1 - T ) r 

^ j — e — ► 1, as fe — > oo, 

r=0 r ' 



where in the middle expression, we write yo = and = x. Hence we have (|2.4|) . 

Part (6) follows since L N = D\ + • • • + -Djv is in the domain of attraction of a stable law ( [HH 
Corollary 2, XVII.5, p. 578]). □ 

We need some additional notation. In this section (Section |2J) we define the giant nodes as the 
k £ largest nodes, i.e., those nodes with degrees D iN) , . . . , D (N _ ks+1) , where k £ is some function of e, 
to be chosen below. We define 

A e ,N = B £jN n C £jJ v n D £ N , (2-5) 

where 

(i) B £jN is the event that the stubs of node 1 and node 2 are attached exclusively to stubs of 
giant nodes. 

(ii) C £jN is the event that any two giant nodes are attached to each other; and 

(iii) D £ ^ N is defined as 

D E , N = {D 1 < b n>£ nD 2 < b DjS }, 
where b Dj£ = min{k : 1 - F(k) < e/8}, so that b Di£ = £ -i/(t-i)(i+o(i)) _ 

The sets B e N and C £)JV depend on the integer k £ , which we will take to be large for e small, and will 
be defined now. The choice of the index k £ is rather technical, and depends on the distributional 
limits defined in Lemma l2.ll Since L N /u N converges in distribution to the random variable rj, with 
support (0, oo), we can find a^ )£ , such that 

F(L N < a V:£ u N ) < e/36, VJV. (2.6) 

This follows since convergence in distribution implies tightness of the sequence L N /u N (jU p. 9]), 
so that we can find a closed subinterval [a, b] C (0, oo), with 

F{L N /u N £ [a,b]) > 1 -e, ViV. 

The definition of b^ k is rather involved. It depends on e, the quantile b D)£ , the value a^ %e defined 
above and the value of r G (1,2) and reads 



i 



2„ \ 2-r 
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where the peculiar integer 2304 is the product of 8 2 and 36. Given b^ ke , we take k e equal to the 
minimal value such that 

P(& £ > hj2) < e/72. (2.8) 

This we can do due to (|2.4() . We have now defined the constants that we will use in the proof, and 
we will next show that the probability of A £ N is at most e: 



Lemma 2.2 For each e > 0, there exists N £ , such that 



£)JV ) <e, N> N £ . 



2.i 



The proof of this Lemma is rather technical and can be found in the appendix. We will now 
complete the proof of Theorem 11.11 sub iect to Lemma 12.21 

Proof of Theorem 11.11 The proof consist of two parts. The event A £tN , implies the event 
{H N < 3}, so that F(A^ N ) < e induces that {H N < 3} with probability at least 1 — e. 

In the first part we show that P {{H N = 1} n A e N ) < e. In the second part we prove that 

lim P (H N = 2) = lim lim P ({H N = 2} n D £ N \ B £ N ) = p, 

N^oo e— >0 N^oo 

for the events B £>N , D £>N defined above and some < p < 1. Since e is arbitrary positive, the above 
statements yield the content of the theorem. 

We first prove that F({H N = 1} n A £)N ) < e for sufficiently large N. The event {H N = 1} 
occurs iff at least one stub of node 1 connects to a stub of node 2. For j < D\, we denote by 
{[l.j] — > [2]} the event that j-th stub of node 1 connects to a stub of node 2. Then, with F N the 
conditional probability given the degrees D\, D2, ■ ■ ■ , -Djvj 



1) A-e,i 



< E 



< E 



J2M{[^j} - [2]}, A, 

3=1 



3=1 



1 



< - < £, 



N - 1 



(2.10) 



for large enough N, since L N > N . 

max{P(^ N ) 1 P(^ Jf )} < P^^) < e 



We prove now that lim P (H N = 2) = p, for some < p < 1. Since by definition for any e > 0, 

N— >oo 



we have that 



\F(H N = 2) - P {{H N = 2}HD £iN \ B £tN ) 



< 



< 



F(H, 



+ 



{H N = 2) - F{{H N = 2}nD £;N n B £tl 



E,N 



F{B £>N ) -1-e 
uniformly in N, for N sufficiently large. If we show that 



lim P ({H N = 2} n D £ . N I B £ , N ) = p T>£ , 



then there exists a double limit 



p T = lim lim P ({H N = 2} n D £ , N \ B £)N ) 



lim P (H N = 2) 

N— >oo 
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Moreover, if we can bound p Tt£ from and 1, uniformly in e, for e small enough, then we also obtain 
that < p T < 1. 

First, to prove the existence of 



lim ¥{{H N = 2}nD £ , N \B £ 

N— >oo 



lim E (Pjv ({H N = 2}HD £N \ B SjN )) = p Tj£ , 

N— *oo 



(2.11) 



we will show that Pjy ({H N = 2} n D £ ^ N \ B £tN ) is a continuous function of the vector 



D 



(N) 



N-k £ +l) 



11. 



This vector, due to l|2.2jl . converges in distribution to (£i, . . . , £fc e , 0). Hence, by the continuous 
mapping theorem A, Theorem 5.1, p. 30], we have the existence of the limit (|2.11|) . We now prove 
the claimed continuity. 

We will show an even stronger statement, namely, that F N {{H N = 2} n D £>N \ B £iN ) is the ratio 
of two polynomials of the components of the vector -Dfc e , where the polynomial in the denominator 
is strictly positive. Indeed, the hopcount between nodes 1 and 2 is 2 iff both nodes are connected 
to the same giant node. For any < i < D±, < j < D% and < k < k £ let Aijk be the event 
that both the i th stub of node 1 and the j th stub of node 2 are connected to the node with the 
(N — k) th largest degree. Then, 



P„ ({H N = 2}nD £ , N \B £!N ) 



Di D 2 ke-1 

U U U A ^ 

i=Xj=l k=0 



B 



where the r.h.s. can be written by inclusion-exclusion formula, as a linear combination of terms 

^ N (An j u ki n • • • n Ai n , jn , kn I b £iN ) . (2.12) 

It is not difficult to see that these probabilities are ratios of polynomials. For example, 



Pat (Ai,j,k I B £ 



D( N -k)(D 



(N-k) 



1) 



(N-ke+l) 

+ --- + D 



(N) 



1) 



(2.13) 



-D(JV-fc) (D(N-k) 



UN 



UN 



1 ' 

U N > 



,D 



(N-k e +l) 
UN 



+ 



ujv '\ u N 



+ 



D 



(AO 



u N 



1 - 

UN ■ 



Similar arguments hold for general terms of the form in (|2.12|) . Hence, ¥ N ({H N = 2} n D £jN \ B £ N ) 
itself can be written as a ratio of two polynomials where the polynomial in the denominator is 
strictly positive. Therefore, the limit (|2.11f) exists. 

We finally bound p £ from and 1 uniformly in e, for any e < 1/2. Since the hopcount between 
nodes 1 and 2 is 2, given -B £)JV , if they are both connected to the node with largest degree, 



¥({H N = 2} n D £jN | B £>N ) > E [F N (A 1MN) | B £:N )} 



and by (|2.13|) we have 

Pr,e 



lim P ({H N = 2}DD £N \B £N )> lim E 

N— >oo ' ' N^oo 



^(JV)P(jyr)-l) 



N)- 



■+D 



(N-k e + l 



)-D 2 



E 



6 



> E 



On the other hand, the hopcount between nodes 1 and 2 is at most 3, given B £tN when all stubs of 
the node 1 are connected to the node with largest degree, and all stubs of the node 2 are connected 



9 



to the node with the one but largest degree. Hence, for any e < 1/2 and similarly to ()2.13|) . we 
have 



Pt 



= lim P ({H N = 2}nD £N \ B £:N ) < 1 - lim P ({H N >2}nD, 

TV^oo N— >oo 

< 1 - lim P ({H N > 2} n Di N | B £ ^ 

TV— >oo \ 2 ' ' 

J-J D (N)- 



S,N I B £jN ) 



< 1 - lim E 

TV^oo 



-2? 



j=0 



D. 



(JV)" 



< 1 



E 



6 



& D,l/2 



(JV-fc e +l)- 



6 



D 2 

n 

i=0 



D 



-2T 



■D(at)H hD( A ,_ fce+1) -D 2 



b _D,l/2 



< 1 



|l«2V D ' 1/2 



Both the upper and lower bound are strictly positive, independently of e. Hence, for any e < 1/2, 
the quantity p TtE is bounded away from and 1, where the bounds are independent of e, and thus 

< p = p T < 1. □ 
This completes the proof of Theorem 1 1 . 1 1 sub i ect to Lemma 12.21 



3 Proof of Theorem HH 

In Theorem ll.2| we consider the hopcount in the configuration model with degrees an i.i.d. sequence 
with distribution given by (|1.4|) . where D has distribution F satisfying (|1.2j) . We distinguish two 
cases: (i) 1/r < a < l/(r — 1) and (ii) a > l/(r — 1). 

We first prove part (ii), which states that the limit distribution of H N is a mixed distribution 
with probability mass p on 2 and probability mass 1 — p on 3, for some < p < 1. Part (ii) of 
Theorem 11.21 is almost immediate from Theorem 11.11 As before we denote by D\ , D<2 , . . . , D N the 
i.i.d. sequence without conditioning, then P(U^ 1 {-Dj > N a }), which is the probability that for at 
least one index i £ {1,2, ... , N} the degree Di exceeds N a , is bounded by 

TV 

J^P(A > = NP{D 1 > N a ) = N 1+a{1 - T) L(N) = N- £ , 

2=1 

for some positive e, because a > l/(r — 1). We can therefore couple the i.i.d. sequence -D (Ar) = 
(D[ N , D ( 2 N) , . . . , ) to the sequence D = (D\,D2, ■ ■ ■ , D N ), where the probability of a miscou- 
pling, i.e., a coupling such that £) (JV) ^ D, is at most N~ £ . Therefore, the result of Theorem 11.11 
carries over to case (ii) in Theorem 11.21 

We now turn to case (i) in Theorem ll.2l We must prove that if we condition the degrees to be 
smaller than iV a with 1/r < a < l/(r — 1), then the graph distance between two arbitrary nodes 
is 3 with probability equal to 1 — o(l). We define 

TV 
n=l 

to be the total degree of the conditioned model. The steps in the proof Theorem ll.2l i) are identical 
to those in the proof of Theorem ll.il however, the details differ considerably. We take an arbitrary 
e > and define an event G £:N such that F(G £ N ) < e, and prove subsequently that for large enough 
N, F({H N = 1} n G £iN ) < e, F({H N = 2} n G £tN ) < e, and that F({H N < 3} n G e>N ) > 1 - e. The 
proof that F(G £ N ) < e is quite technical and moved to the appendix. Since e > is arbitrary, the 
above statements imply that 

lim F(H N = 3) = 1. 

TV^oo 

Let e > be fixed. The event G £)N is defined by 

G £ , N = {D[ N) < b D , £ , DW < b D . £ , if > L £ (N)N 1+a ^}, 
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where, as before, b D>e = min{£; : 1 — F(k) < e/8} and L e (N) is some function satisfying for each 
5 > that 

liminf N S L £ (N) > 1. 

iV^oo 

Lemma 3.1 For each e > 0, there exists N £ , such that, for all N > N £ , 



F(G C £ 



< e. 



(3.1) 



As explained above, the proof of this lemma is rather technical and can be found in the appendix. 

Proof of Theorem ITT21 part (i). We start with the proof that ¥({H N = 1} n G EtN ) < e, for 

sufficiently large N. The event {H N = 1} occurs iff at least one stub of node 1 connects to a stub 
of node 2. As before we denote for j < D^' by {[l.j] — * [2]} the event that j-th stub of node 1 
attaches to a stub of node 2. Then 



P({^ = i}nG E , ff )< E 







< E 


E 










< E 


E 







L 



(N) 



-1{Gs,n} 



b 2 n 

— N-l ^ fc ' 



(3.2) 



for large enough N, since Li N ' > N. 

We will now bound P ({H N = 2} n G S)N ). Note that the event {H N = 2} occurs iff there exists 
node k > 3 and two stubs of node k such that the first one connects to a stub of node 1 and the 
second one connects to a stub of node 2. For k > 3 such that Dj?' > 2, i,j < Dj^\ i ^ j, we 
denote by {[k.i] — > [1], [k.j] — > [2]} the event that the i-th stub of node k connects to a stub of 
node 1 and the j-th stub of node k connects to a stub of node 2. Then, 



: :: ({//v 2}- <-;,,;■ E 



< E 



N Dj» 



Y.Y. ¥ » cum -> w> m ^ [ 2 ]> n 

fc=3 i^j 



N D™ 



D 



(IV) 



D (JV) 



EE r (/V) i r(/V) o 



< ATE((DW) 2 



l{G £iJV } 



(3.3) 



(L e (A0iV 1+a ( 2 - r >-3) 2 ' 



Observe that 

E[(£>W) 2 ] 



£ n<JV a(2n - 1)P(£>1 > < iV Q ) 

PCg^iV) En<jv«(2n - l)L(n - l)(n - T 



l-r 



L 4 (iV)iV a(3 - r) , 



for some slowly varying function L4. Substitution of this upper bound in the right-hand side of (|3,3|) 
shows that for N large enough, 



'({H N = 2}nG e , N ) <NU(N)N 



a(3-r) 



(L e (iV)W 1+a(2_T) - 3) 5 



< e, 



because a(3 — r) < 2(1 + q(2 — r)) when q < 1/(1 — r). 

We will complete the proof by showing that, for sufficiently large N, 



H N < 3}nG £jiV ) > l-e. 



(3.4) 
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Let (3 = (1 + a(4 - r))/4. Since 1/r < a < l/(r - 1), we have 

1 + a( 2 2 " T) < fi < a. (3.5) 
In this section we will call a node k, for 1 < k < N, a giant node if its degree satisfies 

< £)W < iy a_ (3 g) 

We will show below that with probability close to 1 at least one of the stubs of node 1 and at least 
one of the stubs of node 2 are connected to stubs of giant nodes, and that any two giant nodes have 
mutual graph distance 1. This implies that the hopcount between the nodes 1 and 2 is at most 3. 

The non-giant nodes, i.e. nodes with degree less than or equal to N@, are called normal nodes 
. First we will show that the total degree of the normal nodes is negligible with respect to . 
The mean degree of a normal node is 

E[D^1{D < N?}] = F ( D ^ n \ D < N< *) 

n=l 

LTV 3 ] 
n=l 

for some slowly varying function L5. Thus, by the Markov inequality, 



e 

£ 3' 



so that, with probability at least 1 — e/3, the fraction of the contribution from normal nodes on 
G £)JV , for sufficiently large N, is at most 



L e (N)N l +< 2 -r) eL e {N)' 

Since (3 < a, for r S (1, 2) the above fraction tends to 0, as N — > 00. Thus the total contribution of 
the normal nodes is negligible with respect to on G £>N . This implies that for sufficiently large 
N, with probability at most 1 — 2e/3, both nodes 1 and 2 are at graph distance 1 from some giant 
node on G £jJv . It remains to show that any two giant nodes have graph distance 1 with probability 
at least 1 — e/3. For this we need an upper bound of Similarly as above we obtain 

E[D[ N) ] =L e (N)N a ^ 2 - T \ 

for some slowly varying function Lq. Hence from the Markov inequality and since 



\ >s Ls ^ N J " QL 6 (N)N^-r) ~ 6- 



Since we have at most N(N — 1) < N 2 pairs of giant nodes, the probability that two of them, say 
gi and 52 > have graph distance greater than 1 is at most (compare (4.5) of |1 lj ^ . 



+P ^ > h 6 (N)N 1+ < 2 - T ^ < N 2 (l NW 



?L 6 (iV)iV 1 + Q: ( 2 -^ 




2/3-(1+c(2-t)) 

Ar2(e -i + o(ir ^^ + £<£, 
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N=10 3 




N=10 4 




N=10 5 




o 

Q. 



0.2 - 



0.0 



2 3 
hopcount 



Figure 1: Empirical probability mass function of the hopcount for r 
for the unconditioned degrees. 



1.8 and N = 10 3 ,10 4 ,10 5 



for large N, since the exponent grows faster than any power of N. Thus, with probability at least 
1 — |, all giant nodes are on graph distance 1, and thus form a complete graph. This completes 
the proof. □ 



4 Simulation and conclusions 

To illustrate Theorem 11.11 and Theorem 11.21 we have simulated our random graph with degree 
distribution D = |"Z7 — t-i ] , where U is uniformly distributed over (0,1) and where for \x] 
is the smallest integer greater than or equal to x. Thus, 

1 - F(k) = F(U~~^ > k) = jfe 1-1 ", k = 1, 2, 3, . . . 

In Figure E we have simulated the graph distance or hopcount with r = 1.8 and the values of 
N = 10 3 , 10 4 , 10 5 . The histogram is in accordance with Theorem ll.il for increasing values of N we 
see that the probability mass is divided over the values H N = 2 and H N = 3, where the probability 
P(i?jV = 2) converges. 

As an illustration of Theorem 11.21 we again take r = 1.8, but now conditioned the degrees to 
be less than N, so that a = 1. Since in this case (r — 1) _1 = |, we expect from Theorem 11.21 case 
(i), that in the limit the hopcount will concentrate on the value H N = 3. This is indeed the case 
as is shown in Figure El 

Our results give convincing asymptotics for the hopcount when the mean degree is infinite, 
using extreme value theory. Some interesting problems remain, such as 

(i) Can one compute the exact value of p T ? 

(ii) What is the limit behavior of the hopcount when we condition the degrees on being less than 
N a , with a < 1/r? 
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I I N=10 3 




1 2 3 4 5 

hopcount 



Figure 2: Empirical probability mass function of the hopcount for r = 1.8 and N = 
10 3 , 10 4 , 10 5 , 10 6 , where the de grees are conditioned to be less than N. 
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A Appendix. 

A.l Proof of Lemma [2.21 

In this section we restate Lemma 12.21 and then give a proof. 
Lemma A. 1.1 For each e > 0, there exists N e such that 

F(A c £tN ) < e, N > N £ . (A.l.l) 

Proof. We define the event E e ^ N by 



E e,N = E D (N-n) < 8feF7 L ^ 

U=o 

n{D (N _ ke+1) > a Che u N } (b) 
n {L N < b Vfi u N } , (c) 



(A.1.2) 



where b D E is the e-quantile of F used in the definition of D e N and where a^ k ,b^ £ > are defined 
by 

p (&. < <*€ J < e/24 and P(r? > b VjE ) < e/24, 

respectively. Observe that a^ kc is a lower quantile of £fc e , whereas b^ ks defined in (|2.7j) and (|2.8|) is 
an upper quantile of £fc e . Furthermore, b Vj£ is an upper quantile of r], whereas a v ^ e defined in (|2.6|) 
is a lower quantile of r/. Since 

A £) jv = B £ N n C £j jv n D e , n i 
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(see (|2.5|) and below the proof of Lemma [2~T1 for the definition of A £tN , B £>N , C £jJV and D £<N ). we 
have 

f(ai n ) < p(b c £<n nD £>N n e £jN ) + nc c £ , N n d £>n n £ £iJV ) + P(A; C iiV ) + P(^ c iJV ), (A.1.3) 

and in order to prove the lemma we should show that each of the four terms on the right-hand side 
is at most e/4. 

Since on D e N , the nodes 1 and 2 each have at most b D £ stubs, the first term satisfies 



n=0 / 



¥(Bi N nD £>N nE £iN ) < 2b 

due to point (a) of £? £)JV . This bounds the first term of l)A.1.3j) . 

We turn to the second term of (|A.1.3|) . Recall that C £ N induces that no stubs of at least two 
giant nodes are attached to one another. Since we have at most iV 2 pairs of giant nodes g\ and <?2 , 
the items (6), (c) of E £)N imply 

F(ci N nD £iN nE £ , N )< e In 2 J] (i- Ljf _%_ 1 ) 



j=0 

< N 2 (l ' 



= n 2 (e- 1 + o(l)) ) 2 /( 2fe "^) Ujv 

= ^ (e -i + o(1) )K) J /wi'(^ < f> 

for large enough iV, since the exponent grows faster than any power of N. 
The third term at the r.h.s. of (|A.1.3|) is at most e/4, because 

P(Dl N ) < 2P(D 1 > b D , £ ) < 2e/8 = e/4. 

It remains to estimate the last term at the r.h.s. of (|A.1.3|) . Clearly, 



^(tD {N - n) >s^L N ) (a) 

\n=0 '° / 

+P {D (N _ ke+1) < a^u N ) (b) 

+P{L N >b v>£ u N ). (c) 



(A.1.4) 



We will consequently show that each term in the above expression is at most e/12. Let a„ i£ and 
b^ > be as in 1)2. 6 jl and (|2.7j) . then we can decompose the first term at the r.h.s. of (|A.1.4|) as 



D (N _ n) > — — L N < P (L N < a^ujv) + P Aiv-n) > g7 — a^^iijv (A.1.5) 

Vn=0 C,e / Vn=0 D ' e / 

< P (L N < a V)£ u N ) + P (D (N _ hs+1) > b^u N ) 



N 

£ 



+P ^ Al{A < h k u N ] > — 

a ri,e u N 



From the Markov inequality. 

/ AT 



v - e \ iVE (Di 1{D 1 < be. u N \) 
£A1{A < > < 1 ' • (A.1.6) 



vi=l 
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Since 1 — F(x) varies regularly with exponent r — 1 we have, by |1UI Theorem 1(b), VIII. 9, p. 281], 

E(D 1 l{D 1 <b^ k u N })= a -Hk))< 2(2 -r)b^u N (I -F(b ik u N )), (A.1.7) 

for large enough N. Due to (|2.1|) . for large enough N, we have also 

N (1 - < 2. (A.1.8) 

Substituting (|A.1.7jl and (IA.1.81) in (IA.1.6B . we obtain 

2iV(2-r)6^^(l-F(^ £Ujv )) 



AT 

e 



^2 Di\{Di < b£ kc U N } > — — a rh£ u N \ < 



D,e I 8bD £ UNa V,£ 



4(2 -r)^ (l 



for large enough A r . Since 1 — F(x) varies regularly with exponent r — 1, 

(l-Ffe,,,,)) _ 

AT^oo (1-F(t%)) V 4fc ^ 

Hence the r.h.s. of (|A.1.9|) is at most 



8b D ,e ^ 

for sufficiently large N, by definition of b^ k ^ in (|2.7|) . We now show that the second term of ()A.1,5I) 
is at most e/36. Since D / 'ujv converges in distribution to we find from (|2.8J) . 

P {D (N . k£+1) > b^u N ) < P > 6^/2) + e/72 < e/36, 

for large enough iV. Similarly, by definition of a„ E , in (|2.6() . we have 

P(L N < a V:£ u N ) < e/36. 

Thus, the term (|A.1.4l) fa) is at most e/12. 

The upper bound for (|A,1.4|) fb). i.e., the bound 

P (L> (JV -fc £ +i) < a ? fce u iv) < e / 12 , 

is an easy consequence of the distributional convergence of D iN _ ke+1) /u N to ^ e an d the definition 
of a^ fc£ . Similarly, we obtain the upper bound for the term in l|A.1.4[) (c). i.e., 

F {L N > b v>£ u N ) <e/12, 

from the distributional convergence of L N /u N to 77 and the definition of b^ >e . 

Thus we have shown that P(-E? N ) < e/4. This completes the proof of Lemma 12.21 □ 
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A. 2 Proof of Lemma [3.11 

In this section we restate Lemma \3. II and give a proof. 

Lemma A. 2.1 For each e > 0, there exists N e such that for all N > N £ , 

P(G£„)<e. (A.2.1) 

Proof. Clearly, 

F(Gl N ) < 2P (D[ N) > b Dfi ) + P (L^f> < L £ {N)N 1+a V- T ^ . (A.2.2) 
From the definition of b D)£ £ N, and the string of inequalities, 

F(D[ N) > b D>£ ) < F(Dt > b D , £ ) < e/4, 

we obtain that the first term on the right-hand side of (|A.2.2|) is at most e/2. We show now that 
the second term on the right-hand side of (|A.2.2|) is at most e/2. Since 

DjT = E !{ D ?° * fi> 

j<N a 

we obtain, after interchanging the order of the two summations involved, that 

N [N a \-1 N 

^ = EA W = E E^r^il- (A.2-3) 

1=1 j=l 1=1 

We would like to approximate 

N 

HD^ > 3} by iVP (J3f) > j) . 

i=l 

For this we will use (jSJ Theorem 1.7(i), p. 14]) which states that for a binomial random variable 
X with parameters N and p < 1/2, and a such that 

o(l-p)>12 and < < — , (A.2.4) 

we have 



P(\X -Np\ > a) < y^e"3W. (A.2.5) 
Observe that Pn{]) = P (-^1 — j) i s non-increasing in j > 1 and that lim p N (N — 1) = 0. 

iV — >oo 

We will first complete the proof under an extra assumption on {p N (j)}j>i, an d then prove that 
the assumption indeed holds. The extra assumption reads that for some 1 < k < k N < N a , all 
j £ [k, k N ] and N large enough, we have 

(a) p N (j) < 1/2, 

(b) Np N (j) > 288, 

(c) 2Vjj w (j)/432-alogiV-KX), as iV ^ 00, 1 ' ; 

(d) p N {j)>\W(D 1 >j). 

Then, with a = Np N (j)/12, we have that 

(N \ N 

El^ <3}<Np N {j)/2)\ < F(\J2HDf> <j}-Np N (j))\ >a/2)) 

i=l J i=l 

< 12(ATp JV (i))- 1 /2 e -^(i)/432 ) 
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because from (a) and (b) of (|A.2.6|I . we have a(l - p) = (Np N (j)/12)(l - p N (j)) > Np N (j)/24 > 
288/24 = 12 and a/Np = 1/12. Hence, we find 

' k N N \ / k N N 

HiE > j} > ^0')/2} = i -p U^E < 3} * N p-(j)/ 2 } 

yj=K 1=1 J \j=K »=1 

k N / N \ 

^ 1 - E P UE !{ A'"' < J} > ^(i)/2} 



k 



N 



> 1 - E 12(iVPiv(i)) -1/2 exp(-Np N {j)/ 432) >l-k N exp{-Np N (k N ) /432) 

j=K 

> 1 - exp{alog(iV) - Np N (k N )/432} > 1 - e/2, (A.2.7) 
for large A", due to l\A.2.M b)) and JA.2.6f c^. Given (IA.2.7j) and HA.2M (d)). we have 

fcjv N ,j k N Ar k N fcjy- 1 

^ > EE 1 ^" > 3} > y Ep-w 4^ p( ^ j) = t E ^ t ^)< ( a - 2 - 8 ) 

J = K 2=1 J=K j=K— 1 

with probability at least 1 — e/2. 

We will call a function /(A) s/ou>, if for each 5 > 0, /(iV) > A _<5 , for large enough A. Observe 
that for any a > and slow / 

k 

J2f(3)f = h(k)k a+1 , 

3=1 

for some slow function f\. We further assume that we can take k N = Li(A)A Q for some slow 
function L X (N), then the r.h.s. of (IA.2.81) is L £ (N)N 1+a( - 2 - T \ for some slow function L e {N). This 
would imply that the second term of ()A.2.2|> is at most e/2 and thus completes the proof of Lemma 
IA.2.11 subject to l)A.2.6j) and the fact that k N = Li(N)N a for some slow function L\(N). 

We now specify 1 < k < N a , check points l)A.2.6f a))- l|A.2.6f cO) and demonstrate that we can 
take k N = Li(N)N a , for some slow function Li(N). The only restriction on k is (jATHta)). Just 
take re large enough such that p N {n) < 1/2. Then, since j t— > p N {j) is non-increasing, we obtain 
that PnU) < 1/2 for all j > k. Therefore, ()A.2.6f a)) is satisfied. This introduces k and proves 

We next define k N . Choose 

k N = max{P(£>i > k) > 2(1 - F([N a \ - 1))}. 

k 

This definition gives us point (d) of (|A.2.6|) . Indeed for any j < k N , 

" {Di > j) - P (Di > A a ) 



Pjv (j)=P(L>W>j) 



1 - P(-Di > A c 



> P {D x > j) - P (A. > A Q ) > -P (£>i > j) . (A.2.9) 

Before we check the other items of (|A.2.6|) . we prove that k N = L\(N)N a , for some slow function 
Li(N). We argue by contradiction. Suppose there exists 5 > such that k N < A^ 1-5 )" for all A 
sufficiently large. Then, by definition of k N , we have 

P(L>i > k N + 1) < 2P(Di > [N a \) = L(N)N a ^- T \ 

However, when we compute probability on the l.h.s., then we obtain that 

P(Di > k N + 1) > L 3 (A)A-( T - 1 ^ 1 ^ a , 
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for some slow function L%(N). This gives a contradiction. Hence, for any 5 > 0, k N > N a ^ s >, for 
large enough N. This is equivalent to the statement that k N = L\(N)N a for some slow function 
Li(N). 

Since L\(N) is slow, (|A.2.6f c)) follows from (|A.2.6f d)). since 

yP(A > = jL(k N - l)(k N - if- = L 2 {N)N l -^ T - l \ (A.2.10) 

for some slow function Li{N). Since a < l/(r — 1), the right-hand side of (|A.2.10|) increases as 
a positive power of N, and hence exceeds alog(iV), eventually for sufficiently large N. Finally, 
()A.2.6r b)) follows from (|A.2.6f c)). Thus, we have proved that k, k N defined above satisfy the 
assumptions in ()A.2.6|) and that k N = Li(N)N a for some slow function Li(iV). Therefore, the 
second term in (|A.2.2|) is at most e/2, and we conclude that P(Gg N ) < e. □ 
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